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Abstract 



The fully frustrated XY model on a square lattice is studied by means 
of Monte Carlo simulations. A Kosterlitz-Thouless transition is found at 
Tkt ~ 0.446, followed by an ordinary Ising transition at a slightly higher 
temperature, T c ~ 0.452. The non-Ising exponents reported by others, are 
explained as a failure of finite size scaling due to the screening length associ- 
ated with the nearby Kosterlitz-Thouless transition. 

PACS numbers: 64.60.Cn, 75.10.Hk, 75.40.Mg 
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The critical behavior of two-dimensional fully frustrated XY (FFXY) models has been 
the subject of much interest during the last decade. This is seen in the large number of 
papers in the literature, of which a number are very recent. In spite of this, the controversy 
of the nature of the phase transition(s) in these models is by no means settled. 

The models under discussion include the anti-ferromagnetic XY model on a triangular 
lattice ||,||, the square-lattice version of the XY model with one anti-ferromagnetic coupling 
per plaquette [§,|J, and the corresponding Coulomb gas with half-integer charges |5|U. Also 
discussed in this context are the coupled XF-Ising system ]7| and the 19-vertex version of 
the fully frustrated XY model ||, which are believed to be in the same universality class. 
In the present Letter we focus on the square-lattice version of the FFXY model, but since 
the results are expected to have a more general validity, we will repeatedly refer to studies 
of the other above-mentioned models. 

Beside the theoretical questions regarding the universality class, the study of the fully 
frustrated models is largely motivated by their relevance for Josephson junction arrays in 
a magnetic field. Due to this relation, the Hamiltonian of the FFXy model on a square 
lattice is customarily written with the vector potential A^, 

H = -jJ2 C0S (8i ~ 6j + A ij)- 

(ij) 

Here i and j enumerate the lattice sites, 8i is an angle associated with site i, and the 
sum is over nearest neighbors. The frustration is determined by the A^. Full frustration 
corresponds to one half flux quantum per plaquette, which means that / = j- J2 A^ = 1/2, 
where the sum is taken around a plaquette, cf. Eq. (§) below. 

The ground state for this model on a square lattice pj, has plaquettes with clockwise 
and counter-clockwise rotation in a checkerboard pattern. The angular difference between 
nearest neighbors is 4>ij = ®% — 8j + A^ = ±7r/4. This checkerboard pattern gives rise to 
the discrete Z 2 symmetry of the anti-ferromagnetic Ising model, beside the rotational XY 
symmetry. At low temperatures this model therefore has both the topological long-range 
order of the XY model, and the ordinary long-range order of the anti-ferromagnetic Ising 
model. As the temperature is increased, both the XF-like and the Ising-like orders are 
expected to vanish. 

In the first Monte Carlo (MC) study of this model, Teitel and Jayaprakash 0] found 
a steep drop in the helicity modulus, signaling the loss of XY order, accompanied by an 
increase in the specific heat with lattice size, consistent with an Ising transition. Being 
unable to determine the precise critical behavior of the FFXY model, the authors put 
forward two possible scenarios: 

(i) As the Ising temperature, T c , is approached from below, the Ising excitations produce 
a steep drop in the helicity modulus. As this quantity approaches the universal value, 
the KT excitations become important producing a universal jump. This occurs before 
the loss of Ising order, Tkt < T c . 

(ii) As T c is approached from below, the Ising excitations give rise to a jump larger than the 
universal value ||. The loss of Ising and XY order take place at the same temperature, 
^kt = T c . 
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Since then, several investigations have been made with the aim to decide between these 
two possibilities. Whereas some of the earliest MC studies not were decisive a large 

number of recent papers [|^ |8|,[T0| |T^f| on FFXY models, have yielded exponents for the Z 2 



transition which differ from the pure Ising ones, suggesting a new universality class, the 
second possibility above. The evidence is, however, not conclusive, since the finite size 
scalings in these papers not are quite satisfactory. 

In this Letter we present some MC analyses that shed new light on the behavior of the 
FFXY models. We first give evidence for an ordinary KT transition. We then demonstrate 
that the presence of the screening length associated with this transition in the region im- 
mediately above T KT , precludes the use of ordinary finite size scaling, and argue that this 
is the reason for the reported non-Ising exponents. We then determine the Z 2 correlation 
length, present evidence for an Ising temperature T c > Tkt, and demonstrate that our data 
are, indeed, consistent with the pure Ising exponent, v = 1. Our MC data is obtained on a 
bunch of workstations, by means of the ordinary Metropolis algorithm. 

The precise determination of the temperature for a Kosterlitz-Thouless transition is a 
difficult task. This is due both to the absence of spectacular peaks, and a logarithmic 
correction that gives problems with ordinary finite size scaling. A way to cope with these 
difficulties, by extracting the finite-size dependence from the Kosterlitz' renormalization 



group equations, was suggested some years ago |15|. The result may be expressed as a finite 



size scaling relation for the helicity modulus |L6[ valid right at the transition temperature 
Tkt, 

1^ = 1 + \ . (1) 

2T KT 2(lnL + / ) 

Here L is the system size, is the helicity modulus for that size, and Iq is a parameter 
to be determined. The successful application of this relation does, however, require some 
care. In the XY model the amplitude of the spin waves change with both temperature and 
lattice size, which makes it necessary to identify the temperature scale of relevance for the 
vortices, the Coulomb gas temperature, T CG [ PT| , |IB| . In the Coulomb gas this is no problem, 
but to obtain the equivalent of the helicity modulus, one has to include another term in 
the Hamiltonian containing the polarization squared |]r9U2"0|| . The data obtained in this way 
does, indeed, fit very well to Eq. ([]]) plf . 

Figure p] shows the result from this kind of fit for the FFXY model. Note that the scaling 
relation is obeyed only for rather large lattices, L > 32. This is not surprising since Eq. 
(H) is expected to be valid only for low renormalized vortex density. A KT transition at a 
low Coulomb gas temperature |J, as in the model under consideration, has to be monitored 
at larger lengths in order to get data from the region sufficiently close to the critical point. 
The fit gives Tgg = 0.12847(4) corresponding to T KT /J = 0.4460(1). We note that this is 
clearly below T/J ~ 0.454, which is a typical value of the Z 2 transition temperature in the 
literature 

For the study of the Z 2 transition it is customary to define the staggered magnetization 



m, 



(2) 
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where the sum is over all the plaquettes of the system, and m is the vorticity. We define 

the vorticity, in terms of the rotation of the current sin^ = sin(0j — 6j + around a 
plaquette 

m = -^=(sin0 12 + sin 023 + sin0 34 + sin0 4 i). (3) 



The normalization factor in Eq. @ is chosen from the zero-temperature value of m, which 
follows from the angular difference = ±7r/4 in the ground state. 

Several recent MC analyses of the Z 2 transition in the FFXY models make use of the 
expected finite size dependence of various quantities at criticality ||^[Hj|L3j . Such methods 



have generally yielded non-Ising exponents, suggesting a single transition in a new universal- 
ity class. One of several different approaches is to make use of properties of the distribution 



function of M at criticality. This has been done both directly in the FFXY model 0] and 
in the corresponding Coulomb gas || . 

The tacit assumption behind these scaling analyses is, however, that the system size is the 
only relevant length at criticality. The presence of a nearby transition with a corresponding 
characteristic length, A, may well invalidate this assumption. This therefore calls in question 
the attempts to determine the critical exponents of the Z 2 transition in the FFXY models 
by finite size scaling at T c . The condition for a successful application of finite-size scaling 
would be L ^> A, which may imply prohibitively large lattices. 

The effect of this additional length is clearly seen in the dependence of M on the boundary 
conditions (BC). Beside the ordinary periodic BC (PBC) we make use of fluctuating BC 
(FBC) obtained by introducing phase mismatches across the boundaries in the x and 
y directions as a pair of additional dynamical variables. It is with these BC that the XY 
model corresponds to the CG with periodic BC P, |22| , |20|j . Results for the magnetization 



obtained with these two BC, and L = 16, 64, are shown in Fig. [| As the figure shows there 
is a size- dependent temperature region where M is sensitive to the BC. A comparison with 
the helicity modulus (dashed lines), gives at hand that the difference between the two curves 
vanishes when pa 0. This condition implies L ^> A since the vanishing of T means that 
vortex pairs at distance ~ L are free. 

The dependence of M on the BC is presumably a reflection of the dependence of the 
vortex interaction on the BC [20|, an effect that vanishes if T = 0. The conclusion from this 



figure is therefore that the value of M, at fixed T, is uniquely determined by the system 
size, only if L 3> A. 



The implications for ordinary finite size scaling is illustrated with Binder's cumulant |23 

(M 4 ) 



U 



3(M 2 ) 2 ' 



For the usual cases where finite size scaling works, this quantity is size-independent at the 
critical point, Ul = U* . Both the critical temperature and the exponent u, are then obtained 
from the crossing of Ul for different L. Figure |3| shows this kind of analysis for the FFXY 
model with both PBC and FBC, upper and lower points, respectively. A determination of T c 
from data obtained with PBC, would give T c pa 0.454J, in agreement with previous results 



1 1 — 13 1 . About the same critical temperature is obtained from the data obtained with FBC. 



In both cases there is, however, no unique crossing point independent of the lattice sizes. 
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Beside the lack of an unique crossing point, the problem with this kind of analysis is that 
data for the two different BC suggest different values of U*. It is, however, clear that there 
can be only one correct value for U*. This follows since Ul for sufficiently large L will be 
independent of the BC. 

Our conclusions so far for the Z 2 transition are, firstly, that the reported non-Ising 
exponents are artifacts due to the screening length that invalidates the scaling assumption, 
and, secondly, that in order to obtain reliable data for the analysis of the Z 2 transition one 
has to make use of systems with L ^> A, or, equivalently, T « 0. 

For the determination of the Z 2 correlation lengths around T c , we introduce the correla- 
tion function 



g(r) = (-ly^v (m m r ) , (4) 

where the prefactor takes care of the anti-ferromagnetic structure. 

Below T c this function decays to a finite value in the large-r limit, g(r) — > M 2 . The 
decay to this constant is governed by the correlation length 



9{r) 



M 2 ~e- r/f -. (5) 



Figure f| shows g(r) at Tkt- The data is for both PBC and FBC, at lattices of size L = 64, 
128. From the experience with these BC |22| we expect g(r) for an infinite system to lie 



in between the values for the largest lattices. The solid line in Fig. ^ is from Eq. (|5D with 
M 2 = 0.376 and £_ = 6.1, giving strong evidence that the Z 2 transition takes place at a 
higher temperature, T c > Tkt- 

Granted the existence of two separate transitions, one certainly expects the pure Ising 
exponents for the Z 2 transition. To verify that expectation by MC simulations, we turn to 
the temperature-dependence of the correlation length above T c . In that temperature region 
M vanishes, and the correlation function g(r) behaves as 

g(r) ~ e"^. (6) 

Again, we make use of data for lattices and temperatures such that T « 0, which also turns 
out to be a prerequisite for a good fit to Eq. This means that we are only able to obtain 
reliable values for £ at temperatures down to T/J = 0.472. As an additional complication, 
we also have to consider effects from the temperature dependence of the spin waves. 

The effect of the spin waves on the behavior of the FFXY model seems to be overlooked 
so far. In the analysis of the Z 2 transition, the importance of the spin waves stems from 
the fact that the energy associated with a domain wall becomes smaller with larger spin 
wave amplitude. There are two possible contributions to the temperature-dependence: The 
average value of |m| for a single plaquette changes with temperature, and the bare vortex 
interaction may be affected by the spin waves. 

To appreciate the significance of this effect one may well compare with an Ising model, 
Hj = —KJ2(ij) SiSj, with non-singular temperature-dependences in both the coupling con- 
stant and the magnitude of the spins. In the immediate vicinity of T c these temperature 
dependences may be neglected with impunity, but with data in a larger temperature region, 
one has to resort to an effective temperature variable, T 1 = T/(Ks 2 ). To apply this kind of 
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reasoning to the FFXY model we suggest making use of the expression for the bare vortex 
interaction from Ref. [20]. With m defined as in Eq. (§), the energy associated with the 
vorticity at the origin and r becomes 



-Ebarc(r) = — m G(r)m r , 
where J = J (cos0). This suggests that the relevant temperature scale is 

r' = I 

{\m\fj/J 

Figure |5| shows l/£ plotted against T 1 . This is a demonstration that the temperature 
dependence of £ is, indeed, consistent with the pure Ising exponent, v — 1. The analysis 
also suggests a value of the Ising temperature. The figure gives T*j J = 0.4576(13), slightly 
above T^ T / J « 0.440. In ordinary temperatures this corresponds to T c j J = 0.452(1), which 
is consistent both with Tkt as a lower bound, and the crossing points in Fig. [3] as upper 
bounds for T c . 

In conclusion, we have found ample evidence for two distinct transitions in the FFXY 
model on a square lattice. The KT transition is analyzed by finite size scaling of the 
helicity modulus, the previously obtained non-Ising exponents are explained as a failure of 
the scaling assumption, and, with the identification of the relevant temperature scale, the 
data is consistent with the Ising exponent, v — 1. 

The author thanks Professor Petter Minnhagen and Mats Wallin for critical reading of 
the manuscript. Financial support from the Swedish Natural Research Council through 
Contract No. E-EG 10376-302 is gratefully acknowledged. 
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FIGURES 



FIG. 1. Ttt/2T versus lattice size at = 0.12847 together with the finite size scaling func- 
tion, Eq. (|l]). The good fit is strong evidence for an ordinary KT transition. 

FIG. 2. Staggered magnetization, M, with different BC. Data for PBC lie above the corre- 
sponding data for FBC. The dashed lines are the helicity modulus for the same sizes. The data 
differ at temperatures where T^O. 

FIG. 3. Binders cumulant as a function of temperature for the two different BC and several 
lattice sizes. The difference between data from PBC and FBC (solid and dashed lines, respectively) 
casts doubt on the results from this kind of analysis. 

FIG. 4. Correlation function at Tkt- Data is for PBC and FBC, upper and lower points. The 
finite values at large r shows that this temperature is below T c . 

FIG. 5. The dependence of the Z2 correlation length on the effective temperature. With 
^ 1 / u oc T 1 — T/ , the data is consistent with the pure Ising exponent v = 1. 
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